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ABSTRACT 

It has long been suggested that hehum nuclei in the intracluster plasma can sediment in the cluster gravita- 
tional potential well. Some theoretical estimates for the cores of relaxed clusters predict an excess of helium 
abundance by up to a factor of a few over its primordial value. The intracluster helium abundance cannot be 
measured directly. This presents a significant source of uncertainty for cosmological tests based on the X-ray 
derived cluster quantities, such as the gas mass, total mass, and gas mass fraction, all of which depend on the 
assumed helium abundance. We point out that cluster distances derived by combining the Sunyaev-Zeldovich 
(SZ) and X-ray data also depend on the helium abundance. This dependence can be used to measure the abun- 
dance, provided the distance is known independently. For example, if one adopts the WMAP Hq value, then 
the recent Ho measurement by Bonamente and collaborators, derived from SZ data on 38 clusters assuming a 
primordial helium abundance, corresponds to an abundance excess by a factor of L9 + 0.8 within r ~ \ Mpc 
(using only their statistical errors). This shows that interesting accuracy is within reach. We also briefly discuss 
how the SZ and X-ray cluster data can be combined to resolve the helium abundance dependence for the daiz) 
cosmological test. 

Subject headings: cosmic microwave background — distance scale — intergalactic medium — X-rays: galax- 
ies: clusters 



1. INTRODUCTION 

The majority of helium in the universe, predominantly in 
the form of '*He, was produced during the Big Bang. For the 
WMAP value of Qfc/z^ = 0.0223 +0.0007 (Spergel et al. 2007), 
the standard hot big bang nucleosynthesis model predicts a 
primordial fraction of helium in the total baryonic mass den- 
sity of y = 0.2482 + 0.0007 (Walker et al. 1991; Kneller & 
Steigman 2004). Recent spectral measurements in metal-poor 
extragalactic HII regions give a value within 1 % of this theo- 
retical prediction, with similarly small uncertainties (Izotov et 
al. 2007; Peimbert et al. 2007). Thus, the primordial helium 
abundance appears to be known quite accurately. 

Helium abundance in the hot intracluster medium (ICM) 
may differ significantly from the primordial one. First, addi- 
tional helium comes from the stars. The ratio of star mass to 
ICM mass is higher in the cluster centers, so stellar enrich- 
ment will be stronger there. However, the mass of the primor- 
dial helium in the ICM is comparable to the total stellar mass 
in a cluster, so helium enrichment by stars should not be sig- 
nificant (unlike stellar contribution of heavier elements, which 
are present in the ICM in trace amounts). A much greater in- 
crease of helium abundance in the central regions of clusters 
may be caused by sedimentation of heavy nuclei of the ICM 
in the cluster gravitational potential (Fabian & Pringle 1977; 
Rephaeli 1978; Abramopoulos, Chanan, & Ku 1981; Gilfanov 
& Sunyaev 1984; Qin & Wu 2000; Chuzhoy & Nusser 2003; 
Chuzhoy & Loeb 2004; Ettori & Fabian 2006). The consensus 
of the recent works is that if sedimentation is not suppressed, 
then in a hot cluster undisturbed for several gigayears, the 
relative helium abundance can increase by a factor of 2 or 
more within r < 0.2 - 0.3r2()o, and even more at smaller radii 
(Chuzhoy & Loeb 2004). However, sedimentation can be in- 
hibited by several mechanisms, including tangled magnetic 
fields (which should also suppress diffusion and thermal con- 
duction in the ICM, as seems to be observed, e.g., Ettori & 
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Fabian 2000; VikhUnin et al. 2001; Markevitch et al. 2003), 
gas mixing by cluster mergers and turbulence, and the forma- 
tion of a cluster cool core (Ettori & Fabian 2006), because the 
diffusion rate is a strong function of the temperature. 

While mergers and turbulence should inhibit any contem- 
porary sedimentation, they are unlikely to permanently erase 
a large-scale abundance gradient already present by the time 
of the disturbance. The reason is the ICM in relaxed clusters 
is stratified, with low-entropy gas at the bottom of the grav- 
itational well and higher-entropy gas in the outskirts. Such 
a stable gas distribution should restore itself, and any radial 
abundance gradient with it, shortly after a disturbance, pro- 
vided that small-scale ICM mixing during a merger is inef- 
ficient. Indeed, we do observe radially declining iron abun- 
dance profiles in all relaxed clusters with sufficiently detailed 
X-ray data, with the decline traced from the cluster cores to 
at least r ~ 0.5r2oo (e-g-, Fukazawa et al. 1994; Tamura et al. 
2004; Vikhlinin et al. 2005). While the iron abundance gradi- 
ents are probably caused by enrichment rather than sedimen- 
tation, they should be old enough to have survived a merger 
or two, suggesting that it is difficult to erase an abundance 
gradient. Thus, given the uncertainties in the processes that 
inhibit sedimentation in the ICM, it is unclear how significant 
it would be in a typical cluster 

As pointed out by Qin & Wu (2000) and in later works (and 
summarized in 32. II below), helium abundance affects clus- 
ter quantities derived from X-ray observations, such as the 
cluster total mass, gas mass, and gas mass fraction (see also 
Belmont et al. 2005 for an application to the hot gas in the 
center of our Galaxy). The assumed helium abundance also 
affects abundances of heavier elements derived from their X- 
ray emission lines (e.g., Drake 1998; Ettori & Fabian 2006). 
Most of the current X-ray cluster analyses are restricted to 
bright central regions — precisely those regions that may be 
affected by sedimentation. Unfortunately, helium in the ICM 
is fully ionized and not directly observable by spectroscopic 
means. For this reason, its abundance is unknown and has to 
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be adopted from unrelated measurements, e.g., helioseismol- 
ogy (for a review see, e.g., Lodders 2003). The widely used 
X-ray spectral fitting package XSPEC off'ers a choice of abun- 
dance models with the number density of helium relative to 
hydrogen, 

X= , (1) 

rip 

spanning a range between 0.0792 (Lodders 2003) to 0.0977 
(Anders & Grevesse 1989 and some others). For comparison, 
if one takes the abundances of heavier elements to be 0.3 - 
0.5 solar (as in clusters), the CMB-based primordial helium 
abundance (Spergel et al. 2007) corresponds to x = 0.083. 

The unknown cluster helium abundance is a source of un- 
certainty for X-ray cluster-based cosmology studies. In this 
Letter, I propose a way to measure it. 

2. X-RAY DERIVED QUANTITIES 

2 . L Gas mass fraction 

As mentioned above, a helium abundance in the ICM is im- 
plicitly assumed in most X-ray derived cluster quantities, such 
as gas mass, total mass from the hydrostatic equilibrium equa- 
tion, and their ratio /gas. All of these are being used for cos- 
mological tests (see, e.g., Vikhlinin et al. 2003, Henry 1997, 
and Allen et al. 2004 for the three quantities, respectively), 
with projects underway to use them for "precision cosmol- 
ogy". Below, their dependences on the assumed helium abun- 
dance are written in a form relevant for the X-ray analysis 
(that is, fixing the X-ray and other observables). For simplic- 
ity, I assume a uniform abundance over the region involved. 
Current cosmological tests use relatively small central cluster 
regions (r < rgoo or even r < r25oo), for which the effects of 
sedimentation and enrichment can be significant. 

In a fully ionized intracluster plasma, the number of elec- 
trons per proton is 



= 1 + 2x + Xeh ~ 1 



■ 2x, 



(2) 



where x^h represents electrons from elements heavier than he- 
lium; Xeh ~ 0.005 for the intracluster chemical abundances 
(0.3 - 0.5 solar), and we will neglect it. The mean molecular 
weight of the ICM is 



1 + Ax + nik 
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where and jc/, are the mass and number density contribution 
from elements heavier than helium. For clusters, m/, ^ 0.01 
and can be neglected; xi, « Xei, and is certainly negligible as 
well. Thus, for clarity, we will ignore heavy elements below. 

The total mass of a cluster within a certain radius, derived 
under the assumption of hydrostatic equilibrium (e.g., Sarazin 
1988) and (for our illustrative purposes) isothermality, is 
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yU d log r 

Here we used the fact that the logarithmic density gradient 
does not depend on x under our assumption that helium abun- 
dance is spatially uniform; this relation will be more com- 
plicated if one uses a range of radii where this assumption 
does not hold. The electron temperature T^, derived from the 
shape of the X-ray spectrum, is practically independent of the 
helium abundance. 

Let us now consider the cluster gas mass. Provided the clus- 
ter absolute size is known, it is proportional to the plasma 



density. In X-ray analysis, the plasma hydrogen density is 
derived from the normalization of the cluster continuum spec- 
trum (assuming that emission lines from heavy elements are 
detected and properly modeled). The X-ray continuum lumi- 
nosity is 

Lx oc netipeepil + Ax) - n^^eepil + 4x)(l -I- 2x) (5) 

where e^p is bremsstrahlung emissivity for a pure electron- 
proton plasma, and the factor (1 4- 4x) accounts for the addi- 
tional bremsstrahlung on helium nuclei with charge 2. Again, 
heavier elements add very little to the continuum emission. Lx 
is the observable quantity that one obtains directly from spec- 
tral fitting (and a known distance). Fixing it, we obtain the 
dependence of the derived tip on the assumed helium abun- 
dance: 

tip oc [(1 -h4x)(1 +2x)r^'^. 
The gas mass then depends on x as follows: 

1 1/2 



(6) 



Mgas « rip + 4nHe 
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Note that Mgas and Mtot change with x in the opposite di- 
rections, so their ratio, the gas mass fraction, depends on x 
stronger than either of these quantities: 



(1 +Axfl^ 



gas 



(2 + 3x)(l + 2x) 



111- 



(8) 



We do not know the true value of x and have to assume one to 
calculate /gas; Fig. [T] shows the resulting relative error. Note 
that the above /gas(jc) dependence is weaker than that derived 
by Ettori & Fabian (2006); the difference is due to our fixing 
of the observable quantity Lx when deriving Hp in order to 
mimic the X-ray data analysis. 

Fig. [T] shows that an error by a factor of 2 in the assumed 
helium abundance corresponds to a ~ 15% error in /gas. This 
is comparable to the expected difference between the apparent 
/gas at z = 1 for open and flat Q„, = 0.3 cosmologies. Allen et 
al. (2004) derived /gas values within r < r25oo ~ 0.25r2oo for a 
sample of hot clusters, detected such a difference, and used it 




Fig. 1 . — The relative error in the X-ray derived cluster gas mass fraction, 
/gas, as a function of the en'or in the cluster helium abundance, x (see eq.[8). 
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Fig. 2. — The appai'ent Hg value from a combination of SZ and X-ray 
clu ster data, as a function of the error in the cluster helium abundance, x (see 
eq. IllK shown by solid line. For illustration, we take WMAP's Hq = 73 
km s"' Mpc"' as a "true" value (horizontal black dashed Hne; Spergel et al. 
2007). White dashed line and gray en'or band (statistical only, 68%) overlays 
the Hq measurement by Bonamente et al. (2006), who assumed a helium 
abundance from Anders & Grevesse (1989). 
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(Sunyaev & Zeldovich 1972), where the integral is along the 
line of sight, is some characteristic electron density (e.g., 
near the cluster center), and dg is the angular distance to the 
cluster. Here for the rightmost part of the equation, the plasma 
cloud is assumed isothermal and spherically symmetric. The 
latter assumption is important (the line-of-sight distribution of 
the gas density is taken to be the same as that in the plane of 
the sky) but cannot be tested directly, so in practice, a cluster 
sample has to be used to average out any possible ellipticities. 
The surface brightness of the X-ray continuum emission from 
the same line of sight is given by 



{l)np{l)eep{l +4-x)dl oc n\^da 
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The distance can be determined by combining the above equa- 
tions as follows: 



da 
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1 + Ax 
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The quantities y,Sx and T,, are directly measured, but helium 
abundance x has to be assumed. Fig. |2] shows the effect of 
this assumption on the derived dg (or //q oc c/;7'), based on eq. 
(fTTl i. For example, a factor of 2 error in helium abundance 
results in a 10% distance error 



as evidence for Qa > 0. Of course, this cosmological test is 
based on a comparison of /gas at high and low redshifts, so the 
value of helium abundance does not matter as long as it does 
not evolve between those redshifts. Furthermore, the sign of 
the error arising from wrongly assuming a primordial helium 
abundance in the presence of sedimentation is such that the 
derived /gas would be closer to its "pre-sedimentation" value, 
which is what one would ideally want to use for such a test (A. 
Vikhlinin, private communication). A discussion of this error- 
cancellation effect is beyond the scope of this paper (it would 
require a more accurate calculation than that used for eq.|4|i. 
However, it is clear that any tests relying on an even higher 
accuracy of /gas (such as deriving the dark energy equation of 
state, which would need a few percent accuracy on /gas) will 
require the knowledge of the intracluster helium abundance at 
different redshifts. We will return to this in ^4.11 

2.2. Hubble constant 

It has long been suggested (Cavaliere, Danese, & de Zotti 
1979; Silk & White 1978) that combining the X-ray and 
Sunyaev-Zeldovich (SZ) observations of a cluster can be used 
to measure the absolute distance to the cluster This method 
uses the fact that the SZ decrement and the X-ray brightness 
depend on different powers of the intracluster electron density, 
whose value can thus be determined and converted to the dis- 
tance. Below I show how this measurement depends on the 
assumed helium abundance. The underlying reason for this 
dependence is that the SZ effect is caused only by electrons, 
while the X-ray emission is caused by scattering of electrons 
on protons and helium nuclei. The fact that helium sedimen- 
tation would affect the cluster SZ decrement was mentioned 
by Gilfanov & Sunyaev (1984) but, to our knowledge, never 
considered in any SZ data analyses. 

In a nonrelativistic approximation, the SZ signal in the di- 
rection of a cluster is proportional to the comptonization pa- 



3. MEASURING HELIUM ABUNDANCE IN ICM 

The above dependence of the SZ-X-ray distances on the 
helium abundance can be turned around and used to derive 
the cluster helium abundance — provided the distance scale 
is known independently, for example, from Cepheids and su- 
pernovae (Freedman et al. 2001; Sandage et al. 2006) or from 
the CMB fluctuations (Spergel et al. 2007). At present, the 
latter two methods yield similar or smaller uncertainties on 
distances than those from state of the art SZ-X-ray studies 
(e.g., Bonamente et al. 2006). This is because measurement 
errors on the SZ signal are still quite large, while this observ- 
able (y) enters squared in eq. dTTT l. Neither the supernovae 
distances nor those from CMB fluctuations have any signifi- 
cant dependence on primordial helium abundance (Ichikawa 
& Takahashi 2006). 

Fig. |2] shows the value of //q - 77.6;^43 km s"' Mpc"' 
from Bonamente et al. (2006), derived from the SZ and X-ray 
data on 38 clusters at different redshifts. Of their 3 reported 
values, I chose the least model-dependent one, without the 
hydrostatic equilibrium assumption and excluding the central 
cool regions. The above error bars are 68% statistical-only, to 
illustrate an accuracy not quite achieved yet but within imme- 
diate reach (once the instruments have been better calibrated); 
their current systematic uncertainties are twice as big. 

According to M. Bonamente (private communication), they 
assumed the Anders & Grevesse (1989) helium abundance, 
which is a factor of 1 . 17 higher than the primordial value (® . 
Had they used the primordial abundance, their result would be 
//o - 79 km s ' Mpc"^ If, for the sake of argument, we take 
the WMAP value of Ho as "true", and attribute the difference 
between these values to a helium abundance error (Fig.|2]l, we 
conclude that it should be a factor 1 .9 + 0.8 higher than the 
primordial value. Of course, at this accuracy, it is consistent 
with the primordial value, but we can already exclude some 
of the more extreme predictions for helium sedimentation. 
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This and many other current SZ experiments use interfer- 
ometric mapping of the radio brightness. For typical clus- 
ters, this experimental design effectively "subtracts" the sig- 
nal from the shell outside r ~ I Mpc, so the above constraint 
corresponds to this approximate central region. This is about 
the same radius as used for the /gas test by Allen et al. (2004). 

4. SUMMARY AND DISCUSSION 

The usual assumption of a primordial helium abundance 
for the ICM is not necessarily correct, especially in the clus- 
ter central regions where helium may concentrate via several 
mechanisms, such as sedimentation. We cannot measure the 
cluster helium abundance spectroscopically. As pointed out 
by many authors, an incorrect helium abundance assumption 
will result in incorrect cluster gas masses, total masses and 
gas mass fractions derived from the X-ray data. 

In this paper, we point out that cluster distances derived us- 
ing the SZ-X-ray combination also depend on helium abun- 
dance. If one gives up on the original purpose of this method 
and combines it with an independent distance scale estima- 
tor (such as supernovae or CMB fluctuations), one can take 
advantage of this dependence and derive the cluster helium 
abundance. At present, this seems to be the only practical way 
of measuring it. The accuracy may already be interesting — 
if one compares the best SZ-X-ray value for Ho (Bonamente 
et al. 2006) with the CMB value (Spergel et al. 2007), the dif- 
ference corresponds to a helium abundance 1.9 + 0.8 times 
the primordial value within the cluster central 1 Mpc regions 
(68% statistical-only uncertainty). Increasing the number of 
clusters in the sample and improving the SZ data accuracy 
will reduce this uncertainty. 

Chuzhoy & Loeb (2004) proposed another way of con- 
straining the cluster helium abundance — by comparing the 
cluster total masses derived using the hydrostatic equilibrium 
assumption (eq.|4]above) and any other technique independent 
of helium abundance, such as gravitational lensing. However, 
this approach appears less practical at present, because the 
expected apparent mass difference (~ 10% for a factor of 2 
error in helium abundance) is much smaller than persistent 
discrepancies between these mass estimators that are likely 
to be caused by substructure in the dark matter distribution 
and deviations from hydrostatic equilibrium in the ICM (e.g., 
Gavazzi 2005; Meneghetti et al. 2007). In comparison, the 
SZ-X-ray method uses the same object — the ICM — at both 
wavelengths and does not rely on hydrostatic equilibrium. It 
does, however, assume that the ICM is not clumpy (which 



appears to be supported by Chandra imaging), and requires 
accurate mapping of the ICM temperature structure. It also 
needs spherical symmetry, which has always been an issue 
for the SZ method of Ho determination. It can be overcome 
by proper (i.e.. X-ray) selection of a sample of relaxed clusters 
that is also big enough to average out the asymmetries. 

Chuzhoy & Loeb (2004) also pointed out that a higher he- 
lium abundance in the ICM would affect stars that forms out 
of this ICM. Thus, evidence of helium sedimentation may also 
be found in the spectra of the central cluster galaxies. 

4. 1 . The daiz) cosmological test and helium abundance 

The main reason why we want to know the cluster he- 
lium abundances is to remove the related uncertainty from the 
cluster-based cosmological tests, such as the growth of struc- 
ture tests that use cluster mass functions and the daiz) tests 
that use /gas or the SZ-X-ray distances. Obviously, using a 
competing distance estimator to measure helium abundances 
will introduce degeneracies into the resulting cosmological 
constraints, the detailed analysis of which is beyond the scope 
of this paper In principle, the above two distance tests can be 
combined to solve for the dependence on helium abundance, 
because, as seen from eqs. (O and ( fTTT ). the distances derived 
from /gas {da /gas ) and from the SZ-X-ray method depend 
on X differently. One complication is that if one allows for 
helium sedimentation, the basic assumption of the /gas-based 
test, that /gas within a certain central region does not evolve 
with z, may be violated. So in practice, one would have to fit 
together the X-ray and SZ data for a sample of relaxed clus- 
ters spanning a range of z, parameterize and fit any systematic 
change of helium abundance with redshift (hopefully small or 
negligible), and use modeling to deduce "pre-sedimentation" 
/gas values (to within a scaling factor), which would be the 
ones proportional to the universal baryon fraction. To derive 
cosmological parameters other than Hq, the absolute values of 
da (and so the absolute values of x) are not needed, only its 
change with z- However, independent distances for at least a 
few clusters will be required to determine if any helium sedi- 
mentation occurs at all, and if so, to model it. 
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